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Abstract 

Exact relations between the QCD thermal pressure and the trace anomaly are derived. These are used, 
first, to prove the equivalence of the thermodynamic and the hydrodynamic pressure in equilibrium in 
the presence of the trace anomaly, closing a gap in previous arguments. Second, in the temporal axial 
gauge a formula is derived which expresses the thermal pressure in terms of a Dyson-resummed two- 
point function. This overcomes the infrared problems encountered in the conventional perturbation- 
theory approach. 
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CN ■ 1 Introduction 

in 

Conventional methods of calculating the QCD thermal pressure in perturbation theory have encoun- 
tered an apparently insuperable infrared problem W. In previous papers, we have suggested a resum- 
mation technique that may overcome this problem and tested it on scalar field theory^]. Although 
in principle this technique can be applied to QCD, in practice this would be far from simple, since it 
involves introducing a mass for the gluon. Although this mass must disappear at the end of the cal- 
culation, its presence at intermediate stages is a considerable inconvenience. In this paper, therefore, 
we introduce a more suitable technique for performing the resummation that is needed to overcome 
the infrared problem. 

It is common to distinguish two different definitions of pressure^. The first, called the thermodynamic 
pressure, is given in terms of the grand partition function 

z = Y,(i\eM-mm (i.i) 

i 

The summation is over a complete set of physical states i. The Hamiltonian H is an integral over an 
element of the energy-momentum tensor: 

H= f cf^xT 00 ^) (1.2) 



The thermodynamic pressure is defined to be 

P = TV~ l \ogZ (1.3) 
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The other definition of the pressure is in terms of the thermal average of another element of the 
energy-momentum tensor. The hydrodynamic pressure is 

V = (T 33 ) (1.4) 

Although the two pressures P and V are generally expected to be equal for a system in thermal 
equilibrium®, it was remarked some time ago by Landsman and van WeertW that the presence of 
the trace anomaly might cause certain complications. 

In section 2 we derive formulae for the derivatives of the thermodynamic pressure P with respect to 
the bare coupling g and to the temperature: see (2.7) and (2.8). In section 3 we replace g with the 
renormalised coupling g^, and we use the renormalisation group to show that the two formulae for 
the derivatives of the thermodynamic pressure P together imply that P is equal to the hydrodynamic 
pressure V. We show also that the pressure's departure from simple T 4 behaviour is calculable from 
the QCD trace anomaly. In section 4 we show, by working in the temporal gauge, how simple Dyson 
resummation of the thermal gluon propagator removes the infrared problem. In section 5 we derive 
some simple results for the small-coupling limit, which will be the starting point for a perturbation 
expansion. We plan to consider such an expansion in a further paper. 



2 Unrenormalised formulae 



Consider pure-glue QCD for simplicity. At finite temperature there is a choice: one may heat up only 
the physical components of the gauge field® , but for our purposes it will be more convenient to adopt 
the more conventional formalism where also the unphysical components are heated. Then one may 
write® the grand partition function as a path integral: 

/ nr-ifS p 
dAdcdcdB exp i / dx° / cf^x^x) (2.1) 

The Lagrangian is 

£(g,x) = -lF 2 +£ GF + £ GHOS T (2.2) 

We have chosen to implement the gauge fixing by means of an auxiliary field B. The two common 
choices of r are r = 0, with the x° integration running along the imaginary axis, which is the imaginary- 
time formalism, and r = — oo, with the x° integration following the Keldysh contour, which is a version 
of the real-time formalism®. Either formalism may be used for most of our work, though when a 
definite choice has to be made we will choose the second. So far, the fields are unrenormalised. 

We scale each field by some power of g, so obtaining new fields *A, *c, *c, *B. The integration measure 
then acquires a power of g that depends, loosely speaking, on the number of space-time points, so it 
is T-dependent. In order to cancel this, we consider the ratio 

Z(T,0) [ ' 

and make the same field transformation in the denominator as in the numerator, so that the extra 
factor cancels. We then differentiate the logarithm of (2.3) with respect to g, and transform back to 
the original fields. Choosing different powers of g in the definitions of the starred fields and using 
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the fact that the derivative of (2.3) must be independent of what powers we choose, we may obtain 
various identities. For our purposes, the most useful change of field variables is 

A=*A/g c=*c c=*c B =*Bg (2.4) 

We use the definition (1.3) of the thermodynamic the pressure in terms of the grand partition function, 
and the translation invariance property that 

(J ^ dx Q J d n - l KF 2 {x)"j = -i(3V{F 2 (0)} (2.5) 

where V is the volume of the (n — l)-dimensional x-space. We find that 

P(T,g)=P(T,g)-P FREE (T) (2.6) 

satisfies 

§- g P(T,g) = j- g [(F 2 m - <F 2 (0)) FREE ] (2.7) 
Notice the importance of the subtraction term: it avoids a divergence at g = 0. 

We may directly obtain a different derivative of P(T,g). We use the original definition (1.1) of Z, 
together with the expression (1.2) for the Hamiltonian H as an integral over the energy density T 00 : 

§p[PP(T,g)] = -[<T°»> - (T°») FREE ] (2.8) 

where again we have used translational invariance. 

A gauge-invariant form for the gluonic part of the energy-momentum tensor is 

= -F^ p F u p + \g^F 2 (2.9) 

A trace over colour indices is understood. Classically, is traceless, but in quantum theory there 
is an anomaly that changes this. We use dimensional regularisation so that we work intially n = 4 — e 
dimensions. Then 

T\ = ~\eF 2 (2.10) 
Renormalisation leaves behind a nontrivial limit as e — ► 0. 
Introduce the notation 

D M = {F M P pM p) _ (F M pf M )fme (2 n) 

where M is not summed. Then 

^-P(T,g) = ±(D°-(3-e)D 3 ) (2.12a) 
^[(3P(T,g)]=lD° + l(3-e)D 3 (2.126) 

in n = 4 — e dimensions. 
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3 Renormalisation 

In (2.12), the coupling g is unrenormalised, and both D° and D 3 are divergent. We now replace g 
with a renormalised coupling g R . For this, we use dimensional regularisation and the MS scheme, so 
that we define the g R by 

g = n e/2 Z(g R )g R {}i) (3.1a) 

where Z(g R ) is a combination of wave-function and vertex renormalisation factors, and a beta function 
by 

a( \ d 9R 



d/j, 



The beta function has the structure ™ 



P(9r) = ~\^9R + P{9r) 



where for small g R 



Differentiating (3.1a) with respect to ji at fixed unrenormalised coupling g, we obtain 



= 0(g R ) 



d 1 
log Z{g R ) + — 
tdg R g R 



Differentiating it instead at fixed fi, we then have 

dg 



d 1 
— \ogZ(g R ) + — 
9 Log R g R 



dgn 



-dgn 



W{9r) 

The pressure, being a physical quantity, is the same before and after renormalisation: 

P R (T,g R )=P(T,g) 
Of course, P R depends also on \i. Hence (2.7) is equivalent to 



P(9R)-^PR(T,g R ) 



where we have used (2.10). So (2.12a) becomes 

d 



(F 2 (0)} - 



(n > - < 



(3(g R )^-P R (T,g R ) = -\e{D° - (3 - e)^ 3 ) 



(3.16) 

(3.2a) 
(3.26) 

(3.3a) 
(3.36) 
(3.4) 



(3.5a) 



(3.56) 



Because the left-hand side of (3.5) must be finite when e — ► 0, the divergent part of (D° — 3-D 3 ) as 
a function of g R must be of order e _1 exactly. Because the derivative of the pressure in (2.12b) is a 
physical quantity and so must be finite, the divergent parts of D° and D 3 , which are temperature- 
dependent, must be equal and opposite: 



D°(T,g R ) ~ ^9R± + d0(T) gR) + G(( 



D 3 (T,g R )~- 



D(T,g R ) 



+ d 3 (T,g R )+0(e) 



(3.6) 



Again, D° and D 3 depend also on jjl. Hence when e — > 

Um^ ^-P R (T,g R ) = -^—D(T,g R ) = -Lim^ — (3.7a) 
°9r f3(g R ) P{9r) 

or, equivalently, 

P(gR)-^P R (T,g R ) = —D(T, g R ) = (T^> (3.76) 
We have used the fact that < T£ > FREE vanishes in the limit e — > 0. 

We now use the renormalisation group to show that the thermodynamic pressure P is equal to the 
hydrodynamic pressure V defined in (1.4). Introduce the running coupling g(T;g R ) such that 
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dgR 



T—g{T;g R ) = f3{g R )^- g(T; g R ) m g(v, 9r) = 9r (3.8a) 



It is familiar that, to calculate g(T;g R ), one makes use of the identity 

(3(g R )-^g(T;g R ) | t = P(S(T;9r)) (3-86) 

Because P R has dimension T n , we may write 

P R (T,g R ) = T n <t>(g R ,T/ f i) (3.9) 

From renormalisation group theory, 

( f>(g R ,T/^) = ( f>(g(T;g R ),l) (3.10) 

Thus 

T^^r, T/m) = /?(g(T; gfl )) g _ ^ </>(§(T; g R ), 1) (3.11) 
From (3.8b), we may replace f3(g(T; g R ))d/dg(T; g R ) with [3{g R )d/dg R . Hence 

r +1 ^[r^(T, to )] = p {gR )l_p R( T,g R ) (3.12) 

With this equation, together with (2.8) and (3.5a), we obtain 

Pr = -^{[(nc)) - <n(°)>™»] - [<^ »> - <t°») pree ]} 

= ^ T [(^(0))-(T-(0)) FREE ] 

= [(T 33 (0)) - <T 33 (0)) PREE ] (3.13) 

where in the last step we have used the spherical symmetry. From (1.4), this is just the statement 

that the thermodynamic pressure P is equal to the hydrodynamic pressure V . This result is valid in 
spite of the presence of the trace anomaly, so this closes a gap in the conventional proof that explicitly 
ignores itW. 
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This result is valid in spite of the presence of the trace anomaly. However, the anomaly does influence 
the exact form of the pressure. In terms of the quantities introduced in (3.6), in the limit e — > 



p= -i(d° + d 3 -D) (3.14) 

The last term is just 

lD = -iLim^ (n> (3.15) 
We note that (3.12) and (3.5a) give also 

T-^P R (T,g R ) - (4 - e)P R (T, g R ) = [<T£(0)> - (^(0)) FREE ] (3.16) 

so when e — > the departure of the pressure from simple T 4 behaviour is just the result of the QCD 
trace anomaly (The other parts of the energy-momentum tensor, arising from the gauge fixing and 
the ghosts, do not contribute to this trace anomaly®.) 



4 Temporal gauge 



In the gauge A = the expression for D° in terms of the gluon field is particularly simple: 



D° = - 



\(d°A*) 2 ) - ((d°A l ) 2 ) FKEE 



(4.1) 



There are some unanswered questions about the temporal gauge, both at zero temperature^] and more 
particularly in the finite-temperature imaginary-time formalism though the real-time formalism 
may well be free of problems! 11 ]. In the temporal gauge, we can ignore Faddeev- Popov ghosts. 

With the Keldysh contour, 



(4.2) 



D°(0) = J -^^{q«f[A%{q) - A% FHEE (q)} 

thermal matrix propagator^]. This matrix 



where the subscript 12 refers to the element of the 2x2 
has the structure 

A%) = M(g°){(^ - ^4-)A T (q) + ^A L (g) W°) 
where the matrices A7 and Aj, are diagonal and 



e W e y\ q °\ 



M(q°) = 

with n(q°) the Bose distribution (e^ ! — To zeroth order in the coupling 



(4.3) 



(4.4) 



A T (q) 



A T (q) 

A* T (q) 

M?) 

Al(q)_ 



A T (q) 



q 2 + if] 



AL{q) = w 



(4.5) 



where some prescription is needed^] to define the meaning of l/(q ) 2 . As long as a rotationally 
invariant pole prescription is used, the Dyson resummed propagator has the same structure as in 
(4.3), but with 

A T (q) 



qi-n T (q,(3) 

We need not retain the appropriate prescriptions Al in (4.5); they play no part because the self 
energies 11^ and II ^ have zero imaginary parts only at q° = 0, and when we use (4.6) in (4.2) the 
contribution from this point is killed by the (q ) 2 in the numerator of the integrand. When we insert 
(4.2) into (3.7) we obtain 

(3{g R )—P R {T,g R ) 
dgn 

d 4 ~ e Q n , , n.o,. - , n.x f 2 _2_ | 1 1 1 



8e Im / ^^ e tfKW( 1+2n tfH 



q 2 -U T (q,P) q* (q^-U L (q,P) (q ) 2 > 

(4.7) 



The factor of 8 is included because a trace over colour indices was understood in all the equations. As 
in our previous workPl, any infrared divergence problem has been rendered harmless by the Dyson 
resummation: divergences of the self-energies are irrelevant because they now appear in denominators, 
and any zero of the denominators at q = is unimportant because of the powers of q present in d A ~ e q. 

In making use of (4.7) in a semi- or non-perturbative way, which is the goal of subsequent work, we 
have to renormalise accordingly. It is known 

[12] 

that removing a factor Z 2 from the first and third 
terms of the integral (4.7) renders their contribution to the integrand finite. This factor may be found 
by solving the differential equation (3.3a), with the boundary condition that Z(0) = 1. The solution 
may be written in the form 

logZ{9R) + - 2lo <^r)= 2 J d7 T^Mh ( } 



(4.9) 



This integral converges when e — > 0, so the divergent part of Z 2 is 

J 2/3(g R ) 

Z 2 g R e 

So (4.7) finally becomes 

-^P R (T,g R ) = -16 Lim_ Im J -^^-9(q°)(q ) 2 (l + 2n(q )) 

r 2 2 1 1 ~i 

I q 2 - ng(g, J) ~ ? + (g°) 2 - U^(q, 0) ~ W ' ^ 

where II j! and Il£ are convergent thermal self-energies. 



7 



5 Small-coupling limit 



So far, our formulae are exact: we have not used perturbation theory to derive them. We now 
investigate how they behave in the limit of small coupling. We content ourselves with simple results 
that follow without detailed calculation, and leave a more extensive investigation to a future paper. 

Because! 1 ] t 1 ^] p R (T,gn) behaves as — ig^T 4 for small gu, we see from (3.2b) and (3.7b) that 

D = 0{g R ) (5.1) 



We may obtain almost the same information directly from (2.12a), which becomes with (3.6) 

£-Pr(T, g) = U— + d°-?,d 3 -D + 0(e) 
dg 2g L e 



(5.2) 



To lowest order the unrenormalised and renormalised couplings are equal, so to this order the left-hand 
side is again d/dgRpn(T,gji); therefore the singular term on the right-hand side must be at least of 
order g R . Furthermore, in order to give agreement with (3.7b), 



d° -3d 3 = - 



2g R D 



(5.3) 



where the further terms are at least of order g R . 

The pressure may also be calculated directly from (3.14). For this to be compatible with (5.2) in 
lowest order, it must be that d° = d 3 to lowest order. That is 



Pr 



(5.4) 



where the further terms are of order g R at least. 



6 Conclusion 

(3.7) is an exact relation between the QCD thermal pressure and the thermal average of the trace 
anomaly. We have used this to prove nonperturbatively the equivalence of the thermodynamic and 
the hydrodynamic pressure. 

With the plausible assumption that thermal field theory may be formulated consistently in the tem- 
poral gauge, (3.7) can be expressed in terms of a manifestly infrared finite integral (4.10) over the 
full propagator. The absence of a chromomagnetic screening mass poses no particular problem to 
an evaluation of (4.10), so it may provide a suitable starting point for developing new resummation 
techniques. 
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